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^ ' Abstract 

o : 

tH- . The photon and dilepton emission rates from quark gluon plasma 

^^ ' and hot hadronic matter have been evaluated. The in-medium modifì- 

^^ , cations of the particles appearing in the internai loop of the self energy 

diagram are taken into account by using a phenomenological effective 
Lagrangian approach, Brown-Rho and Nambu scaling scenarios. We 
note that the in-medium cffects on the low invariant mass distribution 
of dilepton and transverse momcntum spectra of photon are clearly 
visible. 






Oh! 

^' I. Introduction 



c^ 



k><( \ Numerical simulations of the QCD (Quantum Chromodynamics) equation 

;_( ' of state on the lattice predict that at very high density and/or tempera- 

ture hadronic matter undergoes a phase transition to Quark Gluon Plasma 
(QGP) PI. One expects that ultrarelativistic heavy ion collisions (URHIC) 
might create conditions conducive for the formation and study of QGP. Vari- 
ous model calculations have been performed to look for observable signatures 
of this state of matter. However, among various signatures of QGP, photons 
and dileptons are known to be advantageous, primarily so because, electro- 
magnetic interaction could lead to detectable signal. However, it is weak 
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enough to let the produced particles (real photons and dileptons) escape the 
system without further interaction and thus carrying the inforniation of the 
constituents and their niomentuni distribution in the thernial bath. 

The disadvantage with photons is the substantial background from vari- 
ous processes (thernial and non-thermal) ||2|. Among these, the contribution 
from hard QCD processes is weh understood in the framework of perturba- 
tive QCD and the yield from hadronic decays e. g. vr*^ ^ 77 can be 
accounted for by invariant mass analysis. However, photons from the ther- 
mahzed hadronic gas pose a more difficult task to disentangle. Therefore 
it is very important to estimate photons from hot and dense hadronic gas 
along with the possible modifications of the hadronic properties. 

We organize the paper as follows. In section II we discuss the specific 
reactions which are used to calculate the photon and dilepton emission rate 
from a thermal bath. In section III we discuss the medium modifications 
of hadrons. Section IV is devoted for evolution dynamics. Results are pre- 
sented in section V. 

II. Photon and Dilepton Emission Rates 

The basic aim of URHIC is to distinguish between the following two possi- 
bilities: 

A +A ^QGP— >Mixed Phase^Hadronic Phase 

or 
A + A ^Hadronic Phase 
The former (latter) case where the initial state is formed in QGP (hadronic) 
phase will be called the 'QGP scenario' ('no phase transition scenario'). In 
the present work we will contrast the photon and dilepton spectra originat- 
ing from these two scenarios for SPS and RHIC energies. 

The thermal emission rate of a real photon of energy E and momentum 
p can be expressed in terms of the trace of the retarded photon self energy 
(n^^)at finite temperature [^, Q 

where T is the temperature of the thermal medium. The emission rate of 
thermal dileptons (virtual photon) differs from that of real photon (due to 



difFerent phase space factor) in the following way, 




^"^'imn^^ ' 



^ e,qo/T _ 1 



The photon emission rate due to Compton and annihilation processes 
in QGP was performed in Refs. |||, ^ by applying the Hard Thermal Loop 
(HTL) resummation |0, g]. The rate of hard photon {E > T) emission due 
to these processes is given by [^ 

Recently, the bremsstrahlung contribution to photon emission rate has 
been computed [^ by evaluating the photon self energy in two loop HTL 
approximation. The physical processes arising from two loop contribution 
are the bremsstrahlung of quarks, antiquarks and quark anti-quark annihi- 
lation with scattering in the thermal bath. The rate of photon production 
due to bremsstrahlung process for a two flavor thermal system with E > T 
is given by Q 

and the rate due to q — q annihilation with scattering in the thermal bath 

is given by, 

dR^^P 40 p/^, 

E^^ = ^^ aasET e"^/^ {Jt - Jl) , (5) 

where Jt ~ 4.45 and Jl ~ —4.26. The most important implication of this 
work is that the two loop contribution is of the same order of magnitude 
as those evaluated at one loop |^, ^ due to the larger size of the available 
phase space. The net rate of emission is obtained by adding eqs. ^, § and ^. 
In the hadronic matter (HM) an exhaustive set of hadronic reactions and 
vector meson decays involving vr, p, uj and ry mesons have been considered. 
It is well known Q that the reactions vr/) -^ 1:^,1:7: — > pj , nir — > r/7, 
TTT] — > vr 7 , and the decays p — > vr vr 7 and ùj — > vr 7 are the most important 
channels for photon production from hadronic matter in the energy regime 
of our interest. The rates for these processes could be evaluated from the 
imaginary part of the two loop photon self energy involving various mesons. 



The photon emission rate froni vrp — s- 7r7 via the interniediary ai has also 
been included. 

We have considered quark anti-quark annihilation for the evaluation of 
dilepton emission rate from QGP, which is given by 

dR _ a,q{M) „4^ V Ki{nM/T) 
dM (2^)4 4r " 

with the cross section 

SOvr a^ /7 4m2\ / 2m'^\ ,_, 

For low mass dilepton emission rate we consider the yield from 7r"*"-7r~ — > 



e+e , p — > e+e and co — > e+e (see ||T^, 11| for details). 

III. Medium EfFects 

To study the medium effects on the transverse momentum distribution of 
photons and invariant mass distribution of dileptons from URHIC we need 
two more ingredients. Firstly, we require the variation of masses and decay 
widths with temperature, because the invariant matrix element for pho- 
ton (real and virtual) production suffer in-medium modifications through 
the temperature dependent masses and widths of the participants. As the 
hadronic masses and decay widths enter directly in the count rates of elec- 
tromagnetically interacting particles, the finite temperature effects in the 
cross sections, particularly in the hadronic matter are very important in 
URHIC. 

We consider the following scenarios for the in-medium vector meson 
mass variations in the present work. 
(I) Walecka Model [|I|] 

We have evaluated the in-medium mass of p and co mesons and decay widths 
within the framework of Walecka model. The renormalization procedure of 
Refs. [|^, 14 1 has been used to render the vacuum self energy of vector 



mesons finite. The details of the calculations at non-zero T can be found in 
our previous works [|10| , 11, 15, ^, [l^l and we do not reproduce them here. 



However, from the results of those calculations the variation of nucleon, rho 
and omega masses and the decay width(rp) of rho with temperature at zero 
baryon density can be parametrized as (see also ||l8|) 



m^r/mAT = 1 - 0.0264x^-9^ m*/mp = l-0.l2mx^-^^ (8) 




m* /m,^ = 1 - 0.0438x^-°^ T*/Tp = 1 + 0.664x^ - 0.625x^ (9) 

where asterisk indicates effective mass/width in the medium, x = T/Tc and 
Te = 0.16 GeV. According to our calculation the nucleon, rho and omega 
masses decrease differently. 

(II) Brown-Rho and Nambu Scaling |l|, H] 

In the previous it is mentioned that both the photon and dilepton emission 
rates are proportional to imaginary part of the current-current correlation 
functions (spectral function). We have parameterized the hadronic spectral 
functions with a conspicuous resonance plus a continuum. These spectral 
functions (in vacuum) for the isovector and isoscalar channel have been con- 
strained from experimental data on e~^e~ — > hadrons. At finite temperature 
we have parameterized the vector meson mass (pole) and the continuum 
threshold (wq) of the spectral function as a function of temperature in the 
following way (see Ref. |]1^] for details): 

^ = ^ = (i-^|, (10) 

where A is a sort of dynamical criticai exponent and V stands for vector 
mesons {p and uj). Since the numerical value of A is not known, we take two 
typical cases: A = 1/6 (BR scaling) and 1/2 (Nambu scaling) pO|ì . 

IV. Evolution Dynamics 

The observed photon spectrum originating from an expanding QGP or 
hadronic matter is obtained by convoluting the static (fixed temperature) 
rate, as given by Eqs. (|^) and (^), with expansion dynamics. Therefore, 
the second ingredient required for our calculations is the description of the 
system undergoing rapid expansion from its initial formation stage to the 
final freeze-out stage. In this work we use Bjorken-like |^1[ hydrodynamical 
model for the isentropic expansion of the matter in (1 + 1) dimension. For 
the QGP sector we use a simple bag model equation of state (EOS) with 
two flavor degrees of freedom. The temperature in the QGP phase evolves 
according to Bjorken scaling law T^ t = Tfri. In the hadronic phase we have 
to be more careful about the presence of heavier particles and their change 
in masses due to finite temperature effects. The hadronic phase consists of 
vr, p, a;, r] and ai mesons and nucleons. The nucleons and heavier mesons 
may play an important role in the EOS in a scenario where mass of the 



hadrons decreases with temperature. The energy density and pressure for 
such a system of mesons and nucleons is given by, 



(11) 



i=mesons 



and 



^«=. s (w/'''''é^-(^"^'+(w/''^''4^-*^-^' 



t=mesons 



(27r)3y "3^/""-' " ' (27r)3y "S^jv' 

(12) 
where the sum is over ah the mesons under consideration and N stands for 



nucleons and Ei = Jp^ + mf. The entropy density is given by 

SH = '-^^^ ^ 4aMT)T^ = 4^ffeff(m*(T),r)T3 (13) 

where gcs is the effective statistica! degeneracy. 

Thus, we can visuaHze the effect of finite mass of the hadrons through 
an effective degeneracy gcs{m*{T),T) of the hadronic gas. The variation of 
temperature from its initial value Tj to final value Tf (freeze-out tempera- 
ture) with proper time (r) is governed by the conservation of entropy 

ST = Sin (14) 

The initial temperature of the system is obtained by solving the following 
equation self consistently 

diV. _45C(3)^^2._^3 



dy 27r4 



TTRÌ4a,sTfTi (15) 



where dN^/dy is the total pion multiplicity, Ra is the radius of the system, 
Ti is the initial thermalization time and Ocg = (7r^/90)5feflf- The change in 
the expansion dynamics as well as the value of the initial temperature due 
to medium effects enters the calculation of the photon and dilepton emission 
rates through the effective statistical degeneracy. 

V. Results and Discussions 

Having obtained the finite temperature effects on hadronic properties and 
the cooling laws we now integrate the rates obtained in the previous sections 



over the space-time evolution of the colhsion. We must account for the fact 
that the thermal rates are evaluated in the rest frame of the emitting matter 
and hence the momenta of the emitted photons or dileptons are expressed 
in that frame. Accordingly, the integrai over the expanding matter is of the 

forni 

freeze — out 
dN _ f 4 dR{E*,T{x)) 



formation 

where dV stands for invariant phase space elements: d^p/E for photons and 
d^q for dileptons. E* is the energy of the photon or lepton pair in the 
rest frame of the emitting matter and T{x) is the locai temperature. In 
a fixed frame like the laboratory or the centre of mass frame, where the 
4-momentum of the photon or lepton pair is q^ = {E,q) and the emitting 
matter element d^x moves with a velocity u^ = 7(1, i7), the energy in the 
rest frame of the fluid element is given by E* = u^q^. 

As discussed earlier, ^eff is obtained as a function of T by solving Eq.([T^. 
A smaller (larger) value of g^s is obtained in the free (effective) mass sce- 
nario. As a result we get a larger (smaller) initial temperature by solving 
Eq. (|l5| ) in the free (dropping) mass scenario for a given multiplicity. Naively 
we expect that at a given temperature if a meson mass drops its Boltzmann 
factor will be enhanced and more of those mesons will be produced leading 
to more photons [^. However, a larger drop in the hadronic masses re- 
sults in smaller initial temperature, implying that the space time integrated 
spectra crucially depends on these two competitive factors. Therefore, with 
(without) medium effects one integrates an enhanced (depleted) static rate 
over smaller (larger) temperature range for a fixed freeze-out temperature 
{Tf = 130 MeV in the present case). In the present calculation (Fig. (|l|) 
for Pb+Pb collisions at SPS energies, dN/dy = 600) the enhancement in 
the photon emission due to the higher initial temperature in the free mass 
scenario (where static rate is smaller) overwhelms the enhancement of the 
rate due to negative shift in the vector meson masses (where the initial tem- 
perature is smaller). Accordingly, in the case of bare mass (Nambu scaling) 
scenario the photon yield is the highest (lowest). In case of the Walecka 
model, the photon yield lies between the above two limits. 

In the 'QGP scenario' the photon yield with in-medium mass is lower 
than the case where bare masses of hadrons are considered. However, the 
difference is considerably less than the 'no phase transition scenario'. This 
is because, in this case the initial temperature is determined by the quark 
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Figure 1: Total thermal photon yield corresponding to dN/dy = 600 and 
Ti = 1 fm/c. The solid (long-dash) line indicates photon spectra when 
hadronic matter fornied in the initial state at Tj = 195 MeV (Tj = 220 MeV) 
and the medium effects are taken from Nambu scaling (Walecka model). 
The dotted line represents the photon spectra without medium effects with 
Ti = 270 MeV. The solid (dotted) line with solid dots represent the yield 
for the 'QGP scenario' when the hadronic mass variations are taken from 
Nambu scaling (free mass). 
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Figure 2: Total thermal dilepton yield corresponding to dN/dy = 600 and 
Tj = 1 fm/c. The solid (long-dash) line indicates dilepton spectra when 
hadronic matter formed in the initial state at Tj = 195 MeV (Tj = 220 MeV) 
and the medium effects are taken from Nambu scaling (Walecka model) . The 
dotted line represents the spectra without medium effects with Tj = 270 
MeV. The square (solid dots) line with solid dots represent the yield for the 
'QGP scenario' when the hadronic mass variations are taken from Nambu 
scaling (Walecka model). 
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Figure 3: Thermal photon spectra at RHIC energies. 

and gluon degrees of freedom and the only difference between the two is 
due to the difFerent Ufetimes of the mixed phase. In Fig. (|l]), the photon 
spectra from 'QGP scenario' is compared with that from 'no phase transition 
scenario'; the latter overshines the former. 

The space time integrated dilepton spectra for the 'QGP scenario' and 
'no phase transition scenario' with different mass variation are shown in 
Fig. (pi). The shifts in the invariant mass distribution of the spectra due to 
the reduction in the hadronic masses according to different models are dis- 
tinctly visible. Similar to the photon spectra, the dilepton spectra from 'no 
phase transition scenario' dominates over the 'QGP' scenario for invariant 
mass beyond p peak. 

For RHIC energies, we consider the 'QGP' scenario only as a scenario 
of pure hot hadronic initial state within the format of the model used here 
appears to be unrealistic. The thermal photon yield for RHIC {dN/dy = 
1735) is displayed in Fig. (^. The solid line represents the total thermal 
photon yield originating from initial QGP state, mixed phase and the pure 
hadronic phase. The short dash line indicates photons from quark matter 
(QM) (= pure QGP phase + QGP part of the mixed phase) and the long 
dash line represents photons from hadronic matter (HM) (= hadronic part 
of the mixed phase + pure hadronic phase) . In ali these cases the effective 
masses of the hadrons have been taken from Nambu scaling. For pT > 2 
GeV photons from QM overshines those from HM, since most of these high 
Pt photons originate from the high temperature QGP phase. The dotted 
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Figure 4: Thermal dilepton spectra at RHIC energies. 

and the dotdash lines indicate photon yields froni QM and HM respectively 
with bare masses in the hadronic sector. The HM contribution for the bare 
mass is larger than the effective mass (Nambu) scenario because of the larger 
value of the hfe time of the mixed phase in the earher case. It is important 
to note that for pj- > 2 GeV, the difference in the QM and HM contribution 
in the effective mass scenario is more than the bare mass scenario. 

Thermal dilepton yield at RHIC energies for QGP initial state and for 
different mass variation scenarios are shown in Fig. (^). The shape of the 
peak in the dilepton spectra in case of Walecka model is slightly differ- 
ent (broader) from the other cases because of the larger mass separation 
between p and uj mesons in this case. The dilepton yield beyond the vector 
meson peak is larger in the bare mass scenario because of the larger life time 
of the mixed phase. 

We have evaluated both the real photon and dilepton emission rate from 
QGP and hot hadronic matter by taking info accounts in-medium mass 
modifications of vector mesons according to Walecka Model calculations, 
BR and Nambu scaling. We observe that the in-medium effects on both 
photon and dilepton spectra are clearly visible. 
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